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Abstract 

We analyze with the Bayesian method the solar and KaniLAND neutrino data in 
terms of neutrino oscillations. We show that Bayesian credible regions with a flat prior 
in the tan^ 6i2-Am2i plane strongly support the LMA solution, in agreement with the 
usual chi-square analysis. Other reasonable priors are considered in order to test the 
stability of the LMA solution. We show that priors which favor small or large values 
of the mixing angle lead to minor changes of the allowed LMA region, affecting mainly 
its large-tan^ 612 part. 



1 Introduction 



A fraction of the electron neutrinos produced by thermonuclear reactions in the core of the 
Sun disappear during their travel to the Earth. This deficit is the famous "Solar Neutrino 
Problem", discovered almost forty year ago. There have been many attempts to solve this 
puzzle during the years since its discovery. Some of them were based on modifications of 
the solar model in order to have a lower neutrino production, although there was a conflict 
with the energy spectrum provided by the four first-generation experiments Homestake [1], 
Kamiokande [2], GALLEX/GNO [3] and SAGE [4]. Recent experiments (Super-Kamiokande 
[5] and SNO [6-9]) have improved dramatically the precision of solar neutrino data. The 
Sudbury Neutrino Observatory (SNO) experiment [6-9], in which the fluxes of solar i/e, 
T^e + t^fi + and z/g + 0.15(z/^ + i^r) on the Earth have been measured through charged-current 
(CC), neutral-current (NC) and elastic scattering (ES) interactions, has shown that the solar 
neutrino deficit is due to z/g t transitions. The fit of all solar neutrino data favor the so- 
called Large Mixing Angle (LMA) solution, in which the MSW effect [10] is operative inside 
the Sun. The LMA solution has been confirmed by the reactor experiment KamLAND [11], 
which observed Ue —>■ i^x oscillations of reactor antineutrinos. 

Working in the framework of three- neutrino mixing with a small ^13 (see Ref. [12]), we 
analyzed the solar neutrino data with the effective mixing 

z/e = cos 612 vi + sin 6112 V2 , 

Va = - sin O12 vi + cos O12 V2 , (1) 

where Qyi is the mixing angle and z/i, V2 are the massive neutrinos with masses mi, m2. 
The flavor neutrino Va is a linear combination of z/^ and v^-^ which are indistinguishable in 
solar neutrino experiments and in the KamLAND experiment. In fact, it is easy to show 
that Va — cos6'23 z/^ — sin6'23 Vt (see Ref. [13]), with ^23 — 7r/4 from the Super-Kamiokande 
atmospheric neutrino data [14]. 

In the framework of neutrino oscillations, a standard least-squares analysis method is 
usually adopted to analyze the data in the context of Frequentist Statistics [15,16]. These 
analyses have shown that the solar neutrino problem is solved by the LMA solution. How- 
ever, there is a different statistical approach, the Bayesian Probability Theory. The Bayesian 
analysis of solar neutrino data has several advantages over a Frequentist one (see the discus- 
sions in Refs. [17-19]). Bayesian analyses of early solar neutrino data have been presented 
in Refs. [17,20-22]. 

In this paper, we analyze the updated solar neutrino data and the last KamLAND data. 
The solar and KamLAND neutrino experiments are briefly described in section |21 Section |31 
gives a short description of our implementation of the standard least-squares method. In 
section m we present our results for the Bayesian allowed regions corresponding to a flat prior 
in the plane of the oscillation parameters tan^ ^12 and Amli = vrS^ — ml. In section El we 
investigate the stability of the Bayesian allowed regions for other reasonable choices of the 
prior. Conclusions are presented in section IHl 
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2 Solar neutrino experiments 



The first experiment to observe solar neutrinos, was designed and started by R. Davis and 
his collaborators in the late 1960s. This is the famous radiochemical chlorine experiment 
located in the Homcstakc Gold mine in South Dakota [1]. It is made of 615 ton of C2CI4. 
The weak process used for the detection is 

i/e + '^Cl ^ 37Ar + e- , (2) 

which has an energy threshold of 0.814 MeV. Thus, the event rate measured in the Homestake 
experiment is due mainly to solar neutrinos, with small contributions from ^Be and CNO 
neutrinos. Other radiochemical experiments are the gallium experiments GALLEX/ GNO [3], 
located in the Laboratori Nazionah del Gran Sasso in Italy, and SAGE [4], located in an 
underground laboratory at Baksan. The reaction for the detection is 

i/e + ^^Ga ^ ^^Ge + , (3) 

which has an energy threshold of 0.233 McV. Due to such a low energy threshold, the gallium 
experiments can detect the most abundant neutrinos from the initial pp reaction in the Sun. 
The calculated flux in this part of the neutrino spectrum is much less sensitive to the input 
parameters of the SSM, and thus is more rehable, than that of the ^B neutrinos. 

Kamiokande [2] and Super-Kamiokande [5] are water Cherenkov detectors located in the 
Kamioka mine in Japan. In these experiments solar i/g's are detected through the elastic 
scattering process 

i^e + e~ Ue + e" . (4) 

We use the high-statistics data of the Super-Kamiokande experiment, which is a 50 kton 
detector sensitive to ^B solar neutrinos with an energy threshold of about 5 MeV. 

The SNO experiment [6-9] is located in Sudbury, Ontario, Canada. It is a heavy- water 

Cherenkov detector which can make simultaneous measurements of the ^B solar z/g flux and 
the u^^r flux produced by neutrino oscillations through the chargcd-current and neutral- 
current interactions on deutcrons and the elastic scattering on electrons 

z/e + ci p + p +e- (CC), 
I'x + d —>■ p + n + Vx (NCI), 
i/^ + e-^ Vx + e- (ES). 

The first phase of the SNO experiment was carried out from Nov. 2, 1999 to Jan. 15, 
2001 [8] using pure D2O. The data collected during this phase allowed the SNO collaboration 
to publish in 2002 [7] a model-independent evidence of solar — > v,j,,t transitions. In the 

second phase, from July 26, 2001 to Aug. 28, 2003 [9], 2 ton of NaCl were added to the 
D2O target to enhance the detection efficiency of the NC channel. The SNO salt- phase data 
confirmed the results of the D2O phase, with more precision. 

Different from the solar neutrino experiments, KamLAND [11] is an experiment which 
detects the antineutrinos produced by the decay of heavy nuclei in commercial Japanese and 
Korean nuclear reactors. It is a 1000 ton hquid scintillator detector located in the Kamioka 
mine in Japan. KamLAND detects neutrinos though the inverse beta decay process 

i^e + P ^ e"^ n , (5) 
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which has a 1.8 MeV energy threshold. KamLAND is the first experiment which observe! a 
disappearance of reactor z/g's. The results of the KamLAND experiment provided a definitive 
proof in favor of the LMA solution of the solar neutrino problem [23-25]. In our analysis we 
use the 766 ton-years data [11]. 



3 Standard least-squares method 

Following the tradition [15, 16], we first analyzed the data with a standard least-squares 
analysis (often called "x^ analysis"). In our analysis, the least-squares function for each 
type of solar neutrino experiment is 

x^ = J2 {nr - Rf) ^7/ {R7 - Rf) , (6) 

where R^^^ and Rf^ are, respectively, the experimental values and theoretical predictions of 
the observables. The theoretically calculated rates Rf^ depend on the oscillation parameters 
Amli and tan^0i2. In Eq. (jH)), is the inverse of the covariance error matrix built by 
adding in quadrature the statistical and systematic errors, considering mutual correlations. 
We have used the data of the Homestake, GALLEX/GNO, SAGE, Super-Kamiokande and 
SNO solar neutrino experiments. For the global analysis of solar neutrino data, we define 
the solar neutrino least-squares function as 

= ^Cl,Ga + ^SK + ^SNO ' C^) 

where we have separated the least-squares function of the radiochemical experiments, in 
which only the total rate is measured, and those of the Super-Kamiokande and SNO exper- 
iments, which measured also the energy spectra in the ES and CC reactions. The initial 
®B solar neutrino flux is considered as a free parameter to be determined by the fit, mainly 
through the SNO NC data. For the other solar neutrino fluxes, we assume the BP04 Standard 
Solar Model [26] and treat the uncertainties as described in Refs. [17,27-29]. 

In the KamLAND experiment, the survival probability of the reactor electron antineu- 
trinos can be written as 

P{V,^Ve) = 1 - sin^ 2^12 sinM j, (8) 

where L is the reactor-detector distance in meters, Am^ is expressed in eV^ and E,y in 
MeV. The results of the KamLAND experiment are divided into 13 energy bins above the 
threshold. In the analysis of KamLAND data, we use the least-squares function [30] 

i 

where N-'^'^^ and A^/^ are, respectively, the measured and calculated event numbers in each 
energy bin. The calculated event number N^^^ depends on the oscillation parameters Arn^i 
and tan^6'i2. In Eq. ©, <Jsys is the systematic uncertainty and ?7 is a free parameter to be 
determined by the minimization of the least-squares function. 

In our global analysis of the solar and KamLAND neutrino data, the least-squares func- 
tion is 

= X| + X| . (10) 



th 



T^exp 

X/"P) + X/"Pln' ' 



th 



+ 



a 



(9) 
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4 Bayesian analysis 



In the Bayesian approach it is assumed that there is some prior knowledge of the values of the 
parameters to be determined by the data analysis. This prior knowledge must be quantified 
by a function which is called "prior probability distribution function" (see Refs. [17, 20- 
22]). The Bayes' Theorem allows the calculation of the "posterior probability distribution 
function", which quantifies the knowledge of the values of the parameters provided by the 
data viewed in the light of the prior knowledge. 

In the case of the effective two-neutrino mixing in Eq. (^, the parameters to be deter- 
mined with the statistical analysis of the solar and KamLAND neutrino data are tan^^i2 
and Amli. 

The Bayes' Theorem says that the posterior probability distribution function of tan^^i2 
and Amgi is given by 

f, 2n A 2 in p(D| tan^Ou, Am^i,X)p(tan^gi2, Aml^\I) 
p(tan 012, Am2i|D, J) = p(D\I) ' ^ ' 

where p(D| tan^^^i2, Am2i,X) is the likelihood function and p(tan^0i2, Am2i|X) is the prior 
probability function. D represents the data and X represents all the prior general knowledge 
and assumptions on solar and neutrino physics (in the Bayesian Probability Theory all 
probabilities are conditional). The function p(D|X) is the global likelihood, which acts as a 
normalization constant through the constraint J dtan^0i2 dAm2i p(tan^0i2, Am2i|D,X) = 1, 
leading to 

p(D|X) = J dtan^0i2dAm2iP(D|tan^0i2, Am2i,X)p(tan^0i2, Am2i|X) . (12) 

Assuming that the statistical and systematic errors have normal distributions, the likeli- 
hood function is given by 

p(D|tan20i2, Am^-i ,X) = ^= ^= . (13) 

(27r)^s/2^ (2vr)^K/2^^ 

Here = 41 is the number of solar data points and A^k = 13 is the number of KamLAND 
data points. Xg is the solar least-squares function in Eq. (jH)) and Vs is the corresponding 
covariance matrix. is the KamLAND least-squares function in Eq. and Vk is the 
corresponding covariance matrix. 

The prior distribution of the parameters, p(tan^0i2, Am2i|X), is a subjective ingredient 
which quantifies the knowledge on the values of the parameters which is independent from 
the data to be analyzed. In physics, it is generally considered desirable to assume a prior 
distribution which carries as little information as possible on the values of the parameters, 
in order to obtain an unbiased result. In this spirit, we consider a flat prior distribution in 
the tan20i2-Am2, plane, 

p(tan^6'i2, Am2i|X) = const . (14) 

In this case, Eq. (fTTjl becomes 

^. 2n A 2 i-n T\ p(D|tan^gi2, Am^i,X) 

p(tan t^i2, Am2i|D,X) = — . (15) 

dtan^0i2 dAm2i p(D| tan^6'i2, Aml^,X) 
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Figure 1: The 90%, 95%, 99.73% Bayesian credible regions in the tan^^i2-Am2i plane. Left: 
the regions allowed by the recent KamLAND 766 ton- years data (shadowed regions), together 
with the LMA region (empty contours) obtained from the analysis the solar neutrino data. 
Right: the region allowed by a combined analysis of solar and KamLAND neutrino data. 

Using this expression, we obtained, in the tan^0i2-Am2i plane, the credible regions with 
90%, 95% and 99.73% posterior probability shown in Fig. [T] 

In our Bayesian analysis, a credible region is calculated by choosing the smallest region 
over which the integral of the posterior probability distribution function is the given prob- 
ability level. For example, a 90% Bayesian credible region is the smallest region for which 
the posterior probability of the parameters is 90%. 

Let us emphasize that the statistical meanings of Bayesian and Frequentist regions are 
different (see Refs. [17,31-38]). Nevertheless, it is interesting to compare the allowed regions 
obtained with the Bayesian and least-squares methods, especially in view of the fact that 
many scientists (and most human beings) do not know the meaning of Frequentist allowed 
regions and tend to give them a probability content, as if they were Bayesian credible regions. 
Fig. 12 shows a comparison of the Bayesian credible regions with 90%, 95% and 99.73% 
probability with the Frequentist allowed regions at 90%, 95% and 99.73% confidence level 
obtained with a least-squares analysis. One can see that the Bayesian credible regions are 
similar but a little larger than the corresponding allowed regions obtained with the least- 
squares analysis. 

Differences between Bayesian credible regions and the corresponding Frequentist allowed 
regions can be due to a poor statistical quality of the data or to a wrong theoretical model 
which leads to a bad fit of the data. The close similarity of the Bayesian credible regions 
and the corresponding Frequentist allowed regions that we have obtained is a signal of the 
very good statistical quality of the data and a confirmation of the LMA solution of the solar 
neutrino problem. Note, however, that the similarity of the Bayesian and Frequentist results 
depends on the assumption of a flat prior in the Bayesian analysis. In section we will 
change the prior probability distribution function in order to further test the stability of the 
LMA solution. 

It is also useful to calculate the marginal posterior probability distributions of tan^6'i2 
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Figure 2: Comparisons of the Bayesian and least-squares methods. Going from left to right, 
the three figures correspond to 90%, 99% and 99.73% probability (Bayesian) or confidence 
level (least-squares). 



and Am2i, which are given by 

p(tan2^i2|D,X) = y"dAm2^p(tan2^i2,Am^i|D, J), (16) 

p(Am2i|D, J) = y"dtan2^i2p(tan2ei2, Am^^lD, J) . (17) 

These distributions give information on each of the two parameters independently from the 
value of the other. 

The marginal posterior probability distributions corresponding to the credible regions in 
Fig. ^ are shown in Fig. El 

The marginal posterior probability distribution for tan^^i2, which is mainly determined 
by solar data, has only one peak with an approximate skewed- Gaussian shape. The resulting 
value of tan^6'i2 is 

tan^^is = 0.40l|];[{^ (68% probability range) . (18) 

The allowed intervals with 90%, 95% and 99.73% probability, each of which is given by the 
intersections of the distribution function with the corresponding horizontal line in the left 
panel of Fig. El are 

tan^9i2 = [0.33 , 0.48] (90%) , [0.32 , 0.50] (95%) , [0.29 , 0.57] (99.73%) . (19) 

The marginal posterior probability distribution for Am^i has one main peak and two 
small peaks on the sides. These two minor peaks correspond to the regions in the left panel 
of Fig. ^ which are allowed by the KamLAND data and have a partial overlap with the solar 
allowed region. However, from the right panel in Fig. 13 one can see that they are strongly 
disfavored with respect to the main peak, which gives 

Aml-^ = 8.32+°:^^ X 10"^ eV^ (68% probability range) . (20) 

The allowed intervals with 90%, 95% and 99.73% probability are 

Am^i = [7.76 , 8.89] (90%) , [7.73 , 8.91] (95%) , [7.37 , 9.32] (99.73%) x 10" W . 

(21) 
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Figure 3: The marginal Bayesian distributions: left figure is a separated posterior probability 
distribution for tan^6'i2; right figure is the separated posterior probability distributions for 
Amli- Three horizontal lines in each figure show three integrated probabilities: 90%, 95% 
and 99. 73% (from up to down respectively). 

Let us emphasize that the remarkable precision of the determination of tan^6'i2 and Am^i 
in Eqs. ()18j) - ()2Hl represents an impressive success of the solar and KamLAND neutrino 
experiments. 



5 Stability of the LMA region 

The Bayes' Theorem in Eq. flll|) requires a prior probability distribution function. The 
fiat prior that we have adopted in the previous section is often considered as the best non- 
informative choice. In this section, we explore the implications for the Bayesian credible 
regions of different choices for the prior distribution, which may quantify some prior belief 
or knowledge. 

From theoretical considerations, it is know that flavor transitions of solar neutrinos 
can occur over several orders of magnitude of tan^6'i2 and Am2i, through vacuum oscil- 
lations for Am2i < 10^® eV^ and large mixing angles or resonant MSW transitions for 
10"^ eV^ < Aml^ < lO^^eV^ and 10"^ < tan26'i2 < 10. Hence, one may think that a flat 
prior distribution in the log(tan^6'i2)-log(Am2i) plane, 

p(log(tan^ei2), log(Am2i)|J) = const , (22) 

is more appropriate than the flat prior distribution in the tan^^i2-Am|^ plane in Eq. ()14|1 . 
In fact, this prior was adopted in Refs. [17,21,22]. In the tan^^i2-Am2i plane, it roughly 
corresponds to the prior 

p(tan^g,2,Am^,|J)oc ^^^,^^^^^,^ . (23) 

Slight differences between Eq. |221and Eq. EHlis that in the tan^6'i2-Am2]^ plane the credible 
region by Eq. |221is no more iso contours and the probability at the center is shifted from the 
biggest value, and vise versa. 
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Figure 4: The 90%, 95%, 99.73% Bayesian credible regions in the log(tan^6'i2)-log(Am2i) 
plane obtained with the prior in Eq. ()22|). Left: the regions allowed by the KaniLAND 766 
ton-years data (shadowed regions), together with the LMA region (empty contours) obtained 
from the analysis the solar neutrino data. Right: the region allowed by a combined analysis 
of solar and KamLAND neutrino data. 



Fig. m shows the credible regions with 90%, 95% and 99.73% posterior probability ob- 
tained with the prior in Eq. ()22|) . One can see that they are rather similar to the ones in 
Fig. ^ which have been obtained with the prior in Eq. (fT^ . except for a slight preference of 
small values of the parameters due to the prior. However, a comparison of the right panels 
in Figs. HI and m leads to the conclusion that the data are so good that the priors in Eqs. (fT^ 
and produce almost the same result for the credible regions. Hence, in practice the 
choice between the priors in Eqs. ()14|) and ()23p is irrelevant. 

We consider now a prior which is a flat distribution with respect to Amg^ and a normal 
distribution for the parameter tan^6'i2: 

p(tan20i2, Am^, | J) = — ^ e'^*""' Si^-i^?!'^^^ . (24) 

V27rcr 



This case corresponds to a prior belief in tan^^i2 — /i, with an uncertainty cr. The Bayes' 
Theorem in Eq. gives the posterior probability distribution 

U 2. A ^\T^T\ p(D|tan^gi2,Am^„J)e-(^-^^^--^)^/^-^ 
]9(tan 6'i2, Am2i|D, J) = — . (25) 

dtan2^i2 dAm^^ p(D| tan2^i2, Am^^, J) e'^*^^' Q^^-i^? 



Fig. El shows the 90%, 95%, and 99.73% credible regions for = and a = 0.2,0.3,0.5. 
The value fi = corresponds to a prior belief that the mixing angle should be small. One 
can see from the three panels in Fig. El that a small value of a has the effect to exclude the 
large-tan^ Ou part of the credible regions obtained with a flat prior in section |31 On the other 
hand, the low-tan^ 612 part of the credible regions is not affected by the change of prior. This 
is due to the KamLAND data, which exclude small values of tan^6'i2, as shown in the left 
panel of Fig. HI 
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Figure 5: The 90%, 95%, and 99.73% credible regions corresponding to the posterior prob- 
abihty distribution function in Eq. ()25|1 . with fi = 0. The left, middle and right figures have 
been obtained, respectively, with a = 0.2, 0.3, and 0.5. The dotted hues correspond to a fiat 
prior distribution (same as the right figure in Fig. 
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Figure 6: The 90%, 95%, and 99.73% credible regions corresponding to the posterior prob- 
ability distribution function in Eq. ()25|1 . with /i = 1. The left, middle and right figures have 
been obtained, respectively, with a = 0.2, 0.3, and 0.5. The dotted lines correspond to a flat 
prior distribution (same as the right figure in Fig. P). 

Figs. El shows the 90%, 95%, and 99.73% credible regions for /i = 1 and cr = 0.2, 0.3, 0.5. 
The value /i = 1 corresponds to a prior belief in favor of a large mixing angle, close to 
maximal. The leftmost panel in Fig. IHl shows that, with respect to the case of a fiat prior 
considered in section a small value of a leads to an enlargement of the credible regions 
towards large values of tan^ 612, which are allowed by the KamLAND data, as shown in the 
left panel of Fig. ^ The low-tan^ 612 part of the credible regions is affected only mildly by 
the change of prior. 

In conclusion, in this section we have shown that the LMA solution is stable when 
reasonable priors are chosen in place of the fiat one considered in section 01 with a possible 
shrink or enlargement of the large-tan^ 612 part of the credible regions depending on a prior 
belief in favor of a small or large mixing angle. 
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6 Conclusions 



In this paper, we presented the results of a Bayesian analysis of the solar and KamLAND 
neutrino data. We showed that the Bayesian analysis with a flat prior distribution in the 
tan^^i2-Am2i plane leads to an allowed LMA region for — > t'^,T transitions which prac- 
tically coincides with the one obtained with a standard least-squares (x^) analysis. We 
investigated the stability of the LMA allowed region for other reasonable choices of the 
prior. We have shown that the LMA solution is stable against reasonable variations of the 
prior, with a possible shrink or enlargement of the large-tan^ Qyi part of the allowed LMA 
region if one has a prior belief, respectively, in favor of a small or large mixing angle. 
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